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The universality of semiclassical gravity is investigated by considering the behavior of the quanti- 
ties {(j> 2 ) and (T";,}, along with quantum corrections to the effective Newtonian potential in the far 
field limits of static spherically symmetric objects ranging from stars in the weak field Newtonian 
limit to black holes. For scalar fields it is shown that when differences occur they all result from 
the behavior of a single mode with zero frequency and angular momentum and are thus due to a 
combination of infrared and s-wave effects. An intriguing combination of similarities and differences 
between the extreme cases of a Schwarzschild black hole and a star in the weak field Newtonian 
limit is explained. 

There is a long history of calculations of leading order quantum corrections to the effective Newtonian potential 
' using Feynman diagrams for scalar and other loops in the graviton propagator [l|, [3, S 0| ■ Another approach is to 
compute the corrections to the Newtonian potential by solving the semiclassical Einstein equations in the far field 
limit of a static spherically symmetric object when the quantized fields are in the Boulware state Q. For a star in 
the weak field Newtonian limit this has been done for both gravitons [|| and massless scalar fields with arbitrary 
coupling £ to the scalar curvature 0] • It has also been done for massless scalar fields outside of a Schwarzschild black 
^sO \ hole Q. In each case it was found that the leading order correction depends only on the mass of the source. The 
lack of dependence on the structure of the source would seem to provide evidence of some type of universality for 
semiclassical gravity in the far field limit. However, examination of the results for massless scalar fields in 0,(1] shows 
that the corrections to the Newtonian potential only agree in the cases of minimal (£ = 0) and conformal (£ = 1/6) 
coupling. Interestingly these are the only cases we know of where the corrections for a massless scalar field have been 
, computed using the Feynman diagram technique [l], Q . 
CN) ■ In Refs. @,[l[the quantities (</> 2 ) and (T a b) were also computed in the far field limit. The results for a Schwarzschild 
black hole are M M 

\o ' 

T. (T a b)sch = diag [4, 2, -3, -3] + ^ (£ - 1/6) diag [2, -2, 3, 3] (lb) 



> 
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($ 9 )sch = [2 - m ~ 1/6)] (1c) 



where the quantum correction to the effective Newtonian potential is denoted by $> q . The differences between the 
H | results for a star in the weak field Newtonian limit Q and a black hole are 



A <^ 2 > - -4^8* ( 2a ) 
3M 

A < Ta "> = -^5^^- 1 /6)diag[2,-2,3,3] (2b) 

A($ 9 ) = ^«-l/6). (2c) 

As can be see from Eq. ([2]) there is complete agreement between the results only for £ = 0. For £ = 1/6 there is the 
puzzling situation that there is no agreement for (0 2 ) yet there is agreement for (T a b ) and $ g . For all other values of 
£ there is complete disagreement between the two results. 

To investigate the origin of these agreements and disagreements we will use the formalism of Refs. pjj [H| to 
compute the quantities A(</> 2 ) and A(T a {,). The quantum correction to the Newtonian potential <E> g can be obtained 
by solving the semiclassical backreaction equations and thus is closely tied to the value of (T a i,). 
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Because of Birkhoff's theorem it is well known that the geometry outside of any static spherically symmetric 
distribution of matter (which for simplicity we have been referring to here as a star) is the Schwarzschild geometry. 
The general line element for a static spherically symmetric spacetime can be written in the form 



dr 2 

ds 2 = -f(r)dt 2 + — - + r 2 (d9 2 + sin 2 0d(f> 2 ) 
k(r) 



(3) 



For Schwarzschild spacetime / = k = 1 — with M the mass of the black hole. For a static star, / and k will 
consist of different functions of r inside the star but outside of it they will be the same as for a black hole except that 
now M is the mass of the star. 

To compute the values of (</> 2 ) and (T a h) outside of a static star we make use of the method of Ref. [HIHH which 
is an extension of that developed in Refs. [12]. We begin with {4> 2 ). All details are given in Ref. [HI- 

For a static spherically symmetric spacetime the Euclidean metric is obtained using the relation r = it with r the 
Euclidean time. The Euclidean Green's function obeys the equation 



U x G E {x, x') - £R(» G E {x, x') = 



(4) 



with £ the coupling to the scalar curvature R. Formally the real part of the Euclidean Green's function is equal to 
(</> 2 ) in the limit the points come tog ether. Of course Ge is divergent in this limit. If renormalization is done using 
the method of point splitting [13l . Il4j then one uses the DeWitt-Schwinger expansion to compute the renormalization 
counterterms which are then subtracted from Ge before letting the points come together. Writing (4> 2 ) unr en = Ge 
one has 



lim [Re < 



^ 2 )unron - (<^ 2 )ds1 



(5) 



For an arbitrary separation of the points the Euclidean Green's function can be written as 



Ge(x,x') = 



An 2 



du> cos[w(r 



1=0 



l)Pt(co&7)Cu>ePvt(r<)qut(r>) 



(6) 



Here Pi is a Legendre polynomial, cos 7 = cos 6 cos 6' + sm8 s'm9' cos(0 — 0'), and C^e is a normalization constant 
which is determined by the Wronskian condition 



dq^t dp u 



dr dr 

The radial modes p u £ and obey the mode equation 



-q^e 



1 



d 2 S 
dr 2 



2k k df ldk 
V + 2fdr + 2dr~ 



dS_ 
~dr 



r{r - 2M) 



uP_ £(£ + 1) 



5 = 



(7) 



(8) 



where S refers to one of the modes. Note that p^i is fixed by the boundary condition that it be regular either at 
the event horizon of the black hole or the center of the star while q u i> is fixed by the boundary condition that it be 
regular in the limit r — > 00. Since the mode equation has the same form outside of the star as it does outside the 
event horizon of the black hole, q^e is the same for both in these regions. 

When the points are split in the time direction so that the radial and angular points are the same then a high 
frequency version of the WKB approximation for the radial modes can be substituted into Eq. to obtain the 
quantity (0 2 )wKBdiv This quantity has the same ultraviolet divergences as Ge in the limit that the points come 
together. If it is both added and subtracted from Eq. ([5|) then the result is 



^ )ren 



^ ) numeric 



^ } analytic j 



(9) 



with 



/ numeric 



/ analytic 



/WKBdi 



/WKBdi- 
2 >DS 



The quantity (</> 2 ) analytic has been computed analytically [Tlj in an arbitrary static spherically symmetric spacetime 



while 



/ numeric 



must usually be computed numerically. 
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While not truly local, (0 2 ) analytic is effectively local in the coordinate system we are using since it just depends 
on the functions / and k and some of their derivatives evaluated at a given value of the radial coordinate r. Thus 
this piece is exactly the same outside the event horizon of a black hole or outside the surface of a static spherically 
symmetric star. The quantity (</> 2 ) numeric is explicitly nonlocal and depends on the two radial mode functions p u £ and 
q^g. In the region outside the surface of a star it is 



^ ) numeric 



1 f°° 00 / 2 

— 2 / <kj \( 2£ + l)C«dW(r)g W i(r) - -== 



y/r(r - 2M) I 1 - 2M/r 



(10) 



Our goal is to derive Eq. |(3J) and to use the results to understand the reasons for the agreements and disagreements 
discussed above between the results for a black hole and those for a static spherically symmetric star. Thus we want 
to compute the difference 



A 2 )star - (0 2 ; 



blackholc (H) 



in the region exterior to the surface of the star in the case that the masses of the star and black hole are equal. Since 
the spacetime geometry is the same in this region for both the star and black hole, the renormalization counterterms 
must also be the same. Thus A(</> 2 ) is explicitly finite and we expect that its value should be independent of the 
renormalization scheme used provided that the Wald axioms are satisfied |15j . The same applies to the change in the 
stress-energy tensor A(T Q b) discussed below. 

Outside of the star p» and g* must be linear combinations of the Schwarzschild modes psch and qsch- As discussed 
above = qsch- However, since the boundary condition for is fixed at the center of the star 

P* = apsch + Pqsch ■ (12) 

The normalization constants C* and Cs c h are determined from the Wronskian condition (0 using the appropriate 
mode functions. Using (fT2j) it is easy to show that 



«C» = C Sch . (13) 

Then combining Eqs. (fT0)l - (fT3|) one finds that 

A <^> = TS dcj^£+ l)C Sch ^£ql g . (14) 

This formula is exact. The normalization constant Cs c h can be determined from the asymptotic form of the 
Schwarzschild mode functions. At large r these mode functions approach the corresponding modes in flat space 
which are 

Pflat = LU~ e i e (u!r) (15a) 
fflat = u} e+1 k e (ujr) (15b) 

with ig and kg modified spherical Bcssel Functions. Substituting into Eq. |7]) gives 

Csch - - (16) 

for all values of w and £. 

As long as the coefficients in the mode equation © are finite at the surface of the star the mode functions p and q 
and their first derivatives are continuous there. Matching at the surface r — A gives 

a = P*^-P'f^ (17a) 

g Sc hPsch - p Sch qsch 

P = PgchK - PschP* / 17b N 
^SchPSch - Psch^Sch 

where primes denote radial derivatives, and all modes are evaluated at r = A. Substituting Eq. (fT7|) into Eq. (fH)) 
shows that in general there will be corrections due to the presence of the star and that they will depend on the 
structure of the star since does. 



4 



In the large r limit the metric goes to its flat space value and, since the boundary conditions on the modes q w e are 
imposed at r = oo, these modes approach their flat space values. Thus it suffices to use the flat space values for these 
modes in Eq. (fT4|) . They are given in Eq. (|15b|) and are of the form 



9flat = |u/+i e — ' 



C\ C£ 



uor (ujrY 



(18) 



with c\ — 1 and the rest of the constants Cj depending on the value of I. 

Since the integral and sum in Eq. (|14p are finite the order may be interchanged. Then for a given value of I the 
major contribution to the integral comes from small values of to. A power series in inverse powers of r can be generated 
by successively integrating by parts provided that (3^1 /aw is a regular function of to near to — 0. It is then found that 
the leading order contribution comes only from the t — integral and is 

, i 1 'o„ (1Q) 



16irr A aoo 



This is a remarkable simplification and works for all static spherically symmetric stars provided only that the mode 
functions and their first derivatives are continuous across the surface of the star, that (3ujtja u i is at least a C 1 function 
of lo at to = 0, and that the geometry outside the star is the Schwarzschild geometry jig . 

Using the result (fT9|) we can show that there is no correction at leading order if £ = 0. First note that ifw = £ = £ = 
the solution to the mode equation ([5]) which is regular on the event horizon of a black hole or at the center of a 
spherically symmetric star is poo = constant and without loss of generality the constant can be chosen to be 1. It 
follows immediately from Eq. (fT2"|) that 

o-oo = 1 , Aio = . (20) 

Thus A((f> 2 } = if £ = 0. As shown below, this also implies that A(T a ) = in this case, so the leading order behaviors 
of (<fi 2 } and (T a ) do not depend on the type of matter present so long as it is static and spherically symmetric. 
Turning to the calculation of A(T a ), first note that, as for (<fi 2 ) analytic, we can write 

\-T b) — (T ft) numeric b) analytic (^-0 

with (T a b ) analytic the same outside the surface of the star as outside the event horizon of a black hole of the same 
mass. The nonzero components of (T a t,) numer i c are made up of linear combinations of five mode sums called Si ... S5 
in Ref. [Til] . Ignoring the subtraction terms, which are necessary for renormalization, but unnecessary for the analysis 
below, these sums are of the form 

71 ^° £=Q 



with 



2 



S4 = 



qui + VuiV 



ar ar 



S5 = Cujepujequii ■ (23) 



One can find the changes in these sums ASi for a static star by substituting Eq. (|12[) into the above expressions, 
using Eq. (|13p . and taking out the part that would be there in pure Schwarzschild spacetime. This is exactly what 
was done in deriving Eq. pT|) . which is equivalent to AS5. Operationally one simply replaces C u ePui by Csch^jfjW 
in each of the above expressions. 

Then examination of Eq. (4.5) of Ref. [ll| shows that in the far field limit outside of the star 

A(TS) = U + I) AS! + U I) (A5 2 + A A5 3 - ^AS 5 ) (24a) 



A(T r r ) = -lAS 1 + lAS 2 - 1 ^AS 3 + ^AS A + ^AS 5 (24b) 

2r z r 8r z 



A{T"n) = [2(,-^\(AS 1 +AS 2 )+d^AS 3 ~^AS i ~^AS 5 ) . C2-U-) 
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From the definitions of the modes sums AS n it is immediately clear that in general 

AS, = ^ ■ (25) 
dr 

To go further note that since both in the case of a black hole and a star (T a b ) ~ 1/r 5 examination of Eqs. 
shows that these mode sums in the far field limit go like 



AS! = % , AS 2 = % , AS 3 - % , AS 4 = % , AS 5 = % (26) 
Substituting into Eq. ([23]) gives 

C 4 = -3C 5 . (27) 
If the radial derivative of AS4 is computed and Eqs. (jHJ) and (|2~7| are used, then to leading order one finds 

C 2 = -Ci -C 3 + ^C 5 . (28) 



Following the same procedure for AS 2 gives 



C 2 = 5d+3C 3 -^C 5 . (29) 
Eqs. (|28f and (|29|) can be combined to obtain the relations 

C 2 = \c x + - A C 5 , C 3 = -^Ci + C 5 . (30) 

This appears to be as far as one can go with this type of constraint. However, it is possible to go further for the 
specific leading order form of the mode sums. From Eq. (|23|) it is clear that the only difference between ASi and AS5 
is that the integrand of the former has an extra factor of uo 2 . When the same procedure is followed for computing the 
leading order contribution to AS\ that was followed for computing the leading order contribution for A(</> 2 ) = AS5, 
the only difference is a factor of l/2r 2 arising from the difference between dcuexp(—2idr) and duuj 2 exp(—2u!r). 
Thus one finds that C\ = C5/2. Using Eq. (fP9|) to obtain C5 the final result is 

C 5 = 2C X = ~C 2 = 4C 3 = -\c 4 = . (31) 

5 6 lD7raoo 

If Eqs. (|26|) are substituted into Eqs. ((24)) and the relations (|3Tj) are used then it is found that to leading order 

A(T\) = 3000 - 1/6) diag [2, -2, 3, 3] . (32) 
16^000?" 

This immediately shows that for £ = 1/6 the leading order contributions to the components of A(T a b) are zero even 
if the leading order contributions to the mode sums AS1...AS5 are nonzero in this case. Moreover, from Eq. ([2H]) it 
follows that when £ = 0, A(T a b ) vanishes. 

As discussed in the introduction and shown explicitly in Eq. @, there is agreement between the leading order 
behaviors of (<f> 2 ) in the case £ = and of (T a b ) in the cases £ = 0, 1/6. These results provide a mathematical 
explanation for these agreements. 

A potential weak point in the above derivation of the corrections to {<p 2 ) and (T a b ) is the assumption that /^oo/c^oo is 
a C 1 function of uj at lj = 0. We do not have a proof of this, but we have numerically computed the full corrections for 
the case of constant density stars. A previous calculation was done in Ref. [17| "using the approximation developed 
by Page, Brown and Ottewill" [H[. We have obtained the leading order behaviors of the full numerical corrections 
to both ((f) 2 ) and (T a b) in the far field limit We have also computed these corrections by numerically computing 
the quantity /3oo/o!oo which is given in Eq. (JT7J) , and substituting the result into Eqs. (|19p and (|3"2"|) . In every case 
considered there is agreement between these calculations to within the expected numerical error. 

Next we consider the weak field limit. For the metric ([3]) this means that we can write / = 1 + Sf, k = 1 + 5k with 
I Sf I, \Sk\ and their derivatives assumed to be small. If we write = pa a t + Sp, to first order in small quantities the 
mode equation for the case lu — t — becomes 



d 2 Sp + 2 dSp _ d 2 pf iat f25k | 1 dSf + l(Mk\ dp &at + 

dr 2 r dr dr 2 \ r 2 dr 2 dr J dr 



(33) 
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Note that the derivative terms on the right hand side vanish because pa a t = 1. Then the solution to the equation 
which has the correct boundary condition at r = is 



5 P = Z I r ^ f 1 dr ^lR(r 2 ) ■ (34) 
Jo r i Jo 

In the weak field limit the geometry must be nearly flat everywhere. Thus at the surface of the star 

PSch = Pflat + Spsch , gSch = ?flat + <$<?Sch ■ (35) 

For the oj — I — mode psch = 1 so <5psch = [20(. Substituting Eqs. (I34[) and (|55|) into Eq. (TT7I) and then substituting 
the result into Eq. (fT9|) one finds to first order that 



A 

2 



drr 2 R{r) . (36) 



Einstein's equations can be used to write R = 8tt(p — p r — 2pg) with p the density, p r the radial component of 
the pressure and pg its angular component. Thus the leading order correction in the weak field limit depends on 
the equation of state of the matter that makes up the star. In the Newtonian limit the pressure becomes negligible 
compared to the density and one finds 

This shows that the correction now depends only upon the mass of the star. Finally, using this result for AS§ in Eqs. 
(1261) and (EH) one finds 



A < Ta "> = ~T~2~t£ K - V6) diag[2, -2, 3, 3] (38) 

Comparison with Eqs. shows that in the weak field Newtonian limit our results (157)) and (j38[) are in complete 
agreement with those of Ref . [?J • 

One can ask what would happen for other fields. Massless spin 1/2 and spin 1 fields are conformally invariant so 
it is reasonable to expect that the apparent conformal symmetry, which apparently causes the leading order terms 
in the stress-energy to be the same for a star or black hole for the conformally invariant scalar field, will lead to the 
same result for these fields. Since gravitons are not conformally invariant it is not clear what the relationship between 
the stress-energy of gravitons outside of a star and outside of the event horizon of a black hole will be. 

For massive fields the mass provides an infrared cutoff. If we consider massive scalar fields with arbitrary curvature 
couplings the derivation up to Eq. fTJ]) is the same. However, in Eqs. <P~5]) and fT8"l) uj is replaced by = y/cu 2 + m 2 , 
with m the mass of the field. It is easy to show that in this case the corrections are exponentially suppressed. One 
would expect the same behavior for massive fields of higher spin. 

Our results show that, in general, there is a distinction between the far field limit for a strongly gravitating source 
of matter, a weakly gravitating source of matter, and matter which can be treated in the Newtonian limit, even 
if the leading order result for the latter case depends only upon the mass of the star. Nevertheless in most cases 
of physical interest the leading order contributions to the stress-energy for scalar fields and probably for fields of 
higher spin (except possibly gravitons) are independent of the type of matter or its configuration so long as it is 
spherically symmetric. This leads to an apparent universality in the quantum corrections these fields make to the 
effective Newtonian potential. 
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